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Molecular  dynamics  computer  simulations  have  been  used  to  model  monolayer  and 
trilayer  Langmuir-Blodgett  films  formed  by  model  amphiphilic  molecules. 
Several  factors  have  been  identified  that  control  film  structure  and 
stability,  thus  giving  insight  into  optimum  molecular  design.  Head-groups  and 
the  substituent  position  of  hydrocarbon  chains  should  be  designed  to  maximise 
the  stability  of  layer-like  packing.  Molecular  tilt  is  an  important  factor 
since  it  results  in  local  structural  defects  and  roughening  of  interfaces 
which  could  lead  to  formation  of  long  range  inhomogeneities. 

Nximerical  calculations  of  non-linear  optical  response  have  been  performed  for 
model  LB  films  in  which  molecules  are  treated  as  rigid  strings  of  beads  to 
represent  their  size,  shape  and  orientation.  Interactions  between  molecule^s 
are  obtained  as  planewise  dipole  sums  which  prove  to  be  negligible  except 
within  a  layer.  Although  the  local  fields  vary  little  with  tilt  so  that  the 
optical  properties  of  the  film  depend  primarily  on  the  molecular  response 
referred  to  the  film  axes,  the  pattern  of  nonlinear  response  is  sensitive  to 
molecular  tilt. 

In  parallel  with  the  theoretical  work  experimental  studies  have  been  performed 
on  selected  cemdidates  for  practical  LB  films  with  good  optoelectronic 
properties.  Detailed  studies  show  the  stearyl  derivatives  of  DAN,  MAP  and  NPP 
retain  strong  SHG  activity  but  difficulties  were  encountered  in  forming 
satisfactory  films. 
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MOLECULAR  FILMS  FOR  OPTOELECTRONICS 
Statement  of  the  problem  investigated 

This  experimental  and  theoretical  research  program  was  aimed  at  the  optimisation  of  Langmuir- 
Blodgett  films  for  optoelectronic  applications.  The  work  falls  into  three  areas: 

1.  Molecular  dynamics  computer  modelling 

The  aim  here  was  to  understand  how  the  structure  of  multilayer  LB  films  is  controlled  by  the 
intermolecular  forces  between  the  constituent  amphiphilic  molecules,  and  in  this  way  to  provide 
guidelines  for  the  design  of  stable  films  for  optoelectronic  applications.  On  the  basis  of  the 
results  obtained  from  the  simulation  studies  we  have  been  able  to  identify  several  factors  which 
control  the  stability  and  structure  of  LB  films. 


2.  Theoretical  calculations 

The  aim  was  to  determine  how  the  nonlinear  optical  response  of  molecules  that  form  LB  films 
expresses  itself  in  the  nonlineai'  optical  response  of  the  films.  An  understanding  of  this 
problem  can  help  guide  molecular  design  for  optimal  nonlinear  optical  response.  The  problem 
was  addressed  by  performing  theoretical  calculations  on  simplified  model  structures.  A  key 
feature  of  the  method  was  the  use  of  planewise  dipole  lattice  sums  to  treat  the  interactions 
between  layers. 


3.  Synthesis,  deposition  and  characterisation  of  LB  films  formed  by  p-nitroaniline  derivatives. 

Single  crystal  para-nitroaniline  molecules  are  known  to  form  crystals  with  high  non-linear 
optical  response.  Selected  examples  were: 

DAN  -  4-N,N-dimethylamino-3-acetimidonitrobenzene 
MAP  -  2,4-dinitrophenyl-l-alanine  methyl  ester 
NPP  -  (s)-N-(4-nitro-phenyl)-prolinol. 

The  aim  was  to  synthesise  stearyl  derivatives  to  produce  amphiphilic  molecules  SDAN,  SMAP 
and  SNPP  which  can  be  formed  into  Langmuir-Blodgett  films.  Linear  and  non-linear  optical 
properties  of  stable  films  were  then  characterised. 
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Summary  of  major  results 

(Details  of  the  methods  and  techniques  used  and  results  obtained 
are  given  later  in  this  report) 


(a)  Programs  for  molecular  dynamics  simulation  of  multilayer  LB  films  and  for  the 
calculation  of  planewise  dipole  sums  for  layered  structures  have  been  written  and  tested 
against  previous  data. 

(b)  Simplified  but  computationally  tractable  models  have  been  used  which  reflect  the 
essential  molecular  features  of  the  amphiphilic  molecules.  The  head  groups  and  the 
methylene  groups  of  the  tails  have  been  represented  by  single  interaction  sites  with 
appropriately  chosen  paiameters.  This  "first  approximation"  model  has  helped  us  to 
understand  how  various  molecular  factors  control  structure. 

(c)  This  moleculai-  model  has  been  used  to  simulate  and  characterize  monolayers  on  a  flat 
surface  and  this  work  is  now  published  [  M.  Bishop  and  J.  H.  R.  Clarke,  J.  Chem. 
Phys.,  95,  540  (1991)]. 

(d)  We  have  studied  Y-type  nilayer  LB  films  at  200K  and  300K  and  preliminary  reports  of 
this  work  have  been  published.  The  films  were  built  and  equilibrated  layer  by  layer 
with  the  heads  of  the  first  layer  pointing  towards  the  surface.  Each  trilayer  "unit  cell" 
contained  192  chains  with  periodic  boundary  conditions  were  imposed  in  the  directions 
parallel  to  the  surface. 

(e)  Ibe  tails  stay  predominanUy  in  the  all  trans  configuration  at  the  close  packed  densities 
of  LB  films.  As  the  size  of  head  groups  increases  the  available  spare  volume  for  the 
chains  is  filled  by  tilting.  The  mean  tilt  angles  are  in  agreement  with  the  prediction  of 
the  simple  constant  projected  area  model.  The  similarity  between  collective  and  mean 
tilts  is  evidence  for  strong  correlation  between  molecular  orientations  within  a  layer 
although  there  is  little  correlation  in  the  direction  of  tilt  between  adjacent  layers.  This 
may  be  related  to  the  lack  of  orientation  dependent  forces  between  the  head  groups. 

(f)  The  fraction  of  gauche  defects  is  small  but  nevertheless  increases  with  increasing  tilt. 
For  the  hydrophobic  interface  between  chain  ends  this  can  lead  to  chain  interpenetration 
and  interface  roughness.  This  promotes  the  formation  of  defective  regions  in  the  long 
range  lateral  structure. 
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(g)  The  density  profile  across  the  film  reveals  significant  disorder  in  the  head  group 
interface  between  the  second  and  third  layers  which  is  in  broad  agreement  with  recent 
experimental  measurements  using  atomic  force  microscopy  [J.  Gamaes,  D.K. 

Schwatz,  R.  Viswanathan  and  J.A.N.  Zasadrinski,  Nature,  357, 54  (1992)].  This  can 
be  related  to  the  general  metastability  of  these  layer  structures  which,  for  spherical  head 
groups,  results  in  spontaneous  buckling  of  the  interface. 

(h)  Model  film  structures  have  been  devised  that  incoporate  close  packing  of  the  molecules, 
represent  the  elongated  shape  of  the  molecules  and  allow  for  tilt  of  the  molecules  away 
from  the  normal  to  the  film. 

(i)  Calculations  of  planewise  dipole  sums  have  been  performed  for  various  structure. 

Once  the  molecules  become  sufficiently  elongated  (axial  ratio  >5),  fuiiher  elongation 
makes  little  difference.  In  particular  the  interactions  between  molecules  in  different 
layers  prove  to  be  negligible  even  in  markedly  tilted  structure,  which  simplifies  the 
algebraic  theory. 

(j)  Algebraic  expressions  have  been  derived  for  the  local  polarizing  field,  the  macroscopis 
electric  field,  the  linear  susceptibility,  the  relative  permittivity,  the  refractive  indices,  and 
the  quadratic  and  cubic  nonlinear  susceptibilities.  Because  of  the  weak  interactions 
between  planes,  the  results  for  sequences  of  different  layers  are  simple  combinations  of 
those  for  individual  layers. 

(k)  Simplified  values  of  the  molecular  polarisability  and  hyperpolarisability  components 
have  been  used  in  numerical  calculations  of  the  film  properties  in  (k).  The  local  fields 
show  surprisingly  little  dependence  on  molecular  tilt,  so  that  the  susceptibilities  reflect 
mainly  the  effect  of  tilt  on  these  molecular  responses.  Distortions  in  the  layers  have 
little  effect. 

(l)  Tilt  complicates  the  pattern  of  components  of  the  nonlinear  susceptibilities.  In  the 
cubic  susceptibility,  the  "cascading"  contribution  due  to  the  first  hyperpolarizability  P 
can  be  comparable  in  magnitude  to  the  "direct"  one  due  to  the  second  hyperpolarizability 
y,  so  that  the  former  should  not  be  ignored. 

(m)  Synthetic  schemes  and  purification  procedures  were  developed  for  preparing  the  stearyl 
derivatives.  Ail  three  materials  are  crystalline  solids  at  room  temperature. 

(n)  The  crystalline  powders  still  exhibit  SHG  so  that  the  addition  of  stearyl  groups  does  not 
induce  centrosymmetry. 
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(o)  Floating  monolayers  were  characterised  by  p-A  isotherms  and  Langmuir-Blodgett  films 
by  X-ray  diffraction,  UV/visible  spectroscopy  and  optical  measurements 

(p)  SDAN  was  found  to  form  X,  Y  and  Z-type  films  on  silanised  glass  plates  dependent  on 
the  subphase  conditions. 

(q)  Although  SNPP  exhibits  well-defined  p-A  isotherms,  the  transfer  of  films  on  to  solid 
substrates  is  difficult  and  produces  less  reliable  results. 

(r)  For  SMAP  no  conditions  have  been  found  that  give  a  floating  monolayer,  even  in 
mixtures  with  behenic  acid  which  is  strongly  film-forming. 

(s)  For  LB  films  of  SDAN  the  SHG  coefficient  is  about  10%  of  that  of  the  DAN  crystal. 
SNPP  shows  similar  behaviour  relative  to  NPP  whilst  preliminary  measurements  on 
SMAP  show  very  high  conversion  efficiencies  although  there  is  evidence  of 
microcrystaUite  formation. 
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Conclusions 


To  avoid  disordering  effects  in  multilayer  films  the  optoelectronically-active  head 
groups  and  the  position  of  substituent  aikyl  chains  should  be  carefully  designed  to 
maximise  the  natural  stability  for  layer-like  packing.  The  frequent  absence  of  relevant 
experimental  data  emphasises  the  importance  of  atomistic  molecular  mechanics  studies 
for  exploring  this  criterion.  Strong  interactions  can  cause  the  head  groups  to  leave  the 
surface  in  an  attempt  to  form  more  stable  aggregates.  The  effect  is  frustrated  by  the 
chain  interactions  but  can  nevertheless  give  rise  to  considerable  interfacial  disorder. 

The  large  size  of  head  groups  can  lead  to  an  increase  in  chain  disorder  and  consequent 
chain  penetration  effects  which  give  rise  to  roughening  of  the  tail-tail  interface.  These 
effects  could  be  alleviated  by  designing  multiple  chain  amphiphiles. 

Theoretical  calculations  show  that  local  electric  fields  vaiy  little  with  tilt  so  that  the 
optical  properties  of  the  film  depend  mainly  on  the  molecular  response  referred  to  the 
film  axes.  The  pattern  of  nonlinear  reponse  does  however  vaiy  significandy  with 
molecular  tilt.  The  results  form  a  basis  for  treatments  of  more  realistic  models  in  which 
the  response  is  not  uniformaly  distributed  but  is  concentrated  in  the  NLO  chromophore. 
In  combination  with  the  molecular  dynamics  modelling  of  structure,  such  smdies  offer 
a  systematic  approach  to  molecular  design  for  highly  active  LB  films. 

The  preparation  and  purification  of  the  stearyl  derivatives  has  been  performed  for 
several  candidate  molecules.  Although  much  of  the  NLO  activity  is  retained  the 
production  of  stable  multilayer  films  of  high  optical  purity  requires  further  work  on  the 
detailed  design  of  these  molecules  utilising  results  of  the  modelling  and  theoretical 
calculations. 
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Abstract 

Stearyl  derivatives  of  DAN,  MAP  and  NPP  yield  LB  films  with  second-harmonic  generation  (SHG)  coefficients  measured  as  7. 
100  and  20  x  lO'^esu.  Molecular  dynamics  modelling  predicts  molecular  tilts  of  about  20''.  decreasing  slightly  in  successive 
layers.  Model  calculations  show  how  molecular  tilt  affects  film  optical  response. 


1.  Introduction 

This  paper  presents  results  from  a  coordinated  pro¬ 
gramme  to  design,  prepare  and  characterize  LB  films 
with  high  nonlinear  optical  (NLO)  activity.  The  single 
crystal  para-nitroaniline  derivatives  DAN  [1],  MAP  [2] 
and  NPP  [3]  exhibit  encouraging  efficiencies  for  fre¬ 
quency  doubling  of  laser  radiation.  We  have  therefore 
synthesized  further  derivatives  yielding  Langmuir- 
Blodgett  ( LB)  films  with  potentially  high  SHG  efficien¬ 
cies.  and  have  explored  how  molecular  structural 
features  affect  film  structure  and  nonlinear  optical  be¬ 
haviour. 


2.  Materials,  film  deposition  and  characterization 

Suitable  derivatives  must  retain  the  optical  character¬ 
istics  of  the  parent  NLO  molecule  but  yield  ordered, 
homogeneous,  noncentrosymmetric  LB  films.  Incorpo¬ 
ration  of  the  stearyl  moiety  gave  the  compounds 
SDAN,  SMAP,  and  SNPP,  which  achieved  these  aims. 
Even  the  powders  gave  frequency  doubling,  showing 
that  judicious  addition  of  a  long  hydrocarbon  chain  to 
the  active  centre  need  not  render  the  material  NLO 
inactive,  as  happens  with  some  other  derivatives  [4]. 
The  materials  also  demonstrated  good  LB  film  forming 
properties  and  short  cut-off  wavelengths  in  the  visible. 

X-ray  diffraction  spectra  from  multilayer  LB  films  of 
these  materials  on  silanized  glass  substrates  show  a 
large  number  of  dooi  reflections  for  each  material,  indi¬ 
cating  that  they  are  highly  ordered  (Fig.  1).  For  most 
films  the  transfer  ratio  and  large  dow  spacing  exclude  X 
or  Z  type  film  formation,  but  for  SDAN  under  some 
conditions  both  X  and  Z  type  films  appear  to  be 


2  Theta  /  degree 

Fig.  1 .  X-ray  diffraction  spectrum  from  a  Z-type  LB  film  of  SDAN. 

possible.  There  seems  to  be  a  variation  of  molecular 
orientation  for  some  films,  ranging  from  near  vertical  to 
a  tilt  of  about  30°,  depending  on  the  dipping  tempera¬ 
ture  and  pH. 


3.  Optical  properties  and  measurement  of  SHG 

The  unpolarized  UV/visible  absorption  spectra  show 
that  the  cut-off  is  below  450  nm  for  each  material.  For 
SHG  the  incident  light  was  the  1064  nm  wavelength 
from  a  Q-switched  Nd  :  YAG  laser  producing  pulses  of 
2  mJ  energy  and  10  ns  duration,  focused  down  to  about 
2(X)  pm  spot  size  and  p-polarized.  The  p-polarized  har¬ 
monic  signal  was  detected  in  transmission  as  a  function 
of  angle  of  incidence  (Fig.  2)  through  low  wavelength 
pass  filters,  a  532  nm  bandpass  filter  of  1  nm  bandwidth 
and  neutral  density  filters  of  appropriate  optical  den¬ 
sity,  using  a  gated  photomultiplier  tube.  A  nonzero 
signal  is  observed  at  normal  incidence  indicating 
noncentrosymmetry  in  the  plane  of  the  film.  Quartz 
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Fig.  2.  Second-harmonic  signal  as  a  function  of  angle  of  incidence  to 
the  sample  film. 


Maker  fringes  were  used  as  a  calibration  standard  for 
measuring  SHG  d  coefficients.  The  values  obtained  for 
t//10"’  esu  were  7  for  SDAN,  100  for  SMAP  and  20  for 
SNPP. 


4.  Molecular  dynamics  simulations 

Our  simulations  aim  was  to  understand  how  in¬ 
teractions  between  molecules  control  the  structure 
of  multilayer  LB  films.  We  have  extended  previous 
monolayer  studies  [5-7]  to  investigate  Y  type  trilayers 
with  192  chains,  to  assess  how  molecular  orientation  in 
the  first  layer  influences  the  packing  in  subsequent 
layers. 

Despite  their  apparent  chemical  complexity,  the 
molecules  of  interest  can  be  adequately  represented  by  a 
large  single  interaction  site  (the  NLO-active  aromatic 
head  group)  combined  with  a  set  of  17  connected  sites 
(the  CH2  groups  on  the  hydrocarbon  tail).  Lennard- 
Jones  12-6  potentials  were  used  [5,  7]  to  represent  the 
site  interactions.  Values  of  the  Lennard -Jones  parame¬ 
ters  Ific  head-group  mass 

( 100  K,  4.238  nm  and  0.050  kg  mor‘)  were  selected  to 
correspond  roughly  [6]  to  stearic  acid. 

The  chain  was  kept  connected  in  the  dynamics  by 
constraining  the  hydrocarbon  tail  to  0. 1 53  nm  between 
sites  and  (ff^  +  <t,)/2  between  the  head  and  the  first  site. 
The  bond  angle,  torsion  angle  and  surface  potentials 
were  as  in  our  previous  study  [7].  The  first  16  000  time 
steps  (40  ps)  were  discarded  for  equilibration;  the  sub¬ 
sequent  production  run  lasted  40  ps. 

The  mean  tilt  angles  in  the  three  layers  were 
20.62  ±  0.08°,  19.40  ±  0.10°  and  18.42  ±  0.08°  at  a  sur¬ 
face  coverage  of  19.6  per  chain,  very  close  to  that 
obtained  from  our  monolayer  studies  [7], 
20.34°  ± 0.04°.  The  density  G(z)  normal  to  the  r=0 
plane  is  shown  in  Fig.  3.  The  first  layer  is  highly 
ordered,  with  a  sharp  first  peak  arising  from  the  head- 
group  locations.  The  corresponding  peaks  in  the  second 


Fig.  3.  The  density  profile  normal  to  the  r  =  0  plane  for  the  trilayer 
system.  G(:)  is  normalized  so  that  the  area  under  the  curve  for  each 
layer  is  unity.  Distances  are  quoted  as  multiples  of  the  hydrocarbon 
C-C  bond  length. 


and  third  layers  become  increasingly  broader,  suggest¬ 
ing  increasing  roughening  of  the  head-group  layers. 


5.  Theory  of  film  NLO  response 

The  nonlinear  optical  response  of  ordered  films  can 
be  calculated  from  the  molecular  hyperpolarizability  p, 
the  molecular  orientation,  and  the  local  electric  fields  in 
the  layers.  In  turn,  the  local  fields  depend  on  sums  of 
dipole  tensor  interactions  between  a  molecule  in  one 
layer  and  all  the  molecules  in  that  or  another  layer  [8]. 
Calculations  of  these  planewise  dipole  sums  L  have 
been  performed  for  a  series  of  close-packed  model 
structures.  The  molecules  are  represented  as  sphero- 
cylinders  of  diameter  D  and  length  sD  and  are  treated 
as  a  set  of  j  equally  spaced  centres  between  which  the 
interactions  are  averaged.  The  molecules  may  be  ori¬ 
ented  normal  to  the  layers  or  tilted  at  an  angle  6  to  the 
vertical.  Interactions  between  a  molecule  and  those  in 
its  own  plane  dominate  for  all  values  of  s  studied,  with 
interactions  beyond  the  neighbouring  layer  negligible. 
For  5  =  5,  the  in-plane  sums  for  molecules  tilted  to¬ 
wards  the  next-nearest  neighbours  (consistent  with  the 
results  from  the  dynamic  modelling)  vary  as  shown  in 
Table  1,  where  x  is  the  next-nearest-neighbour  direc¬ 
tion.  The  tilt  lowers  the  symmetry  from  hexagonal  to 
monoclinic,  but  olT-diagonal  elements  remain  small. 

For  a  single  layer,  the  local-field  tensor  d  is  given  by 
(1— £.  «)■',  where  a=alsoV,  with  a  the  molecular 
polarizability  and  v  the  molecular  volume.  The  linear 
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TABLE  I.  Lorentz-factor  tensor  components  as  a  function  of 
molecular  tilt  9 


flC”) 

t-v. 

i-v.- 

L,,. 

C-r 

0 

0.7885 

0 

0.7885 

-0.5769 

10 

0.7749 

-0.0052 

0.7851 

-0.5603 

20 

0.7277 

-0.0193 

0.7834 

-0.5112 

30 

0.6299 

-0.0316 

0.8126 

-0.4425 

40 

0.4574 

-0.0018 

0.9670 

-0.4244 

TABLE  2.  Quadratic  susceptibility  components  z.yt/lO  ° 
function  of  molecular  tilt  9 

esu  as  a 

9{°) 

Zvw 

X.V.Vf 

Z,.-.- 

Z,.-.- 

Zi.v.v 

Zzw 

z.-.-.- 

0 

0 

0 

6.0 

5.5 

0 

0 

-17.7 

10 

1.5 

0.4 

9.8 

5.6 

-  4.1 

-1.5 

-15.9 

20 

7.0 

1.6 

I3.I 

5.7 

-10.0 

-3.1 

-13.8 

30 

17.0 

3.7 

14.6 

5.7 

-16.4 

-4.6 

-10.8 

40 

32.5 

6.1 

13.0 

4.7 

-21.6 

-5.4 

-  6.3 

susceptibility  x*"  (whence  the  refractive  indices)  is  given 
by  a  •  d.  We  take  the  molecular  polarizability  to  have 
components  (x,  ot,  jct)  in  molecular  axes,  giving  a  polar* 
izability  per  unti  volume  essentially  independent  of  s  for 
large  s.  With  a/cov  =  1.5  we  find  that  x  and  r  remain 
the  principal  axes  since  L^.  is  so  small.  As  the  molecules 
tilt  towards  x  and  away  from  2,  increases  from  1.07 
to  1.19  at  40°  tilt,  where  n.  decreases  from  1.22  to  1.17; 
tiy  remains  essentially  constant  at  1.07. 

The  quadratic  susceptibility  is  given  by  d"'’  •  b :  dd, 
where  b  =  p/CoV  with  p  the  molecular  hyperpolarizabil¬ 
ity.  We  take  the  molecular  hyperpolarizability  to  be 
essentially  one-dimensional  along  the  molecular  axis 
L,  with  nonzero  components  of  the  form  = 


ifi,  P,  —sp)  in  molecular  axes  A  =  L,  M,  N,  a.  form 
suggested  by  MO  calculations  on  NLO  molecules.  With 
PIcoV  =  1.87  X  10~^®esu  we  find  the  results  for  x,y*  as  a 
function  of  tilt  9  shown  in  Table  2,  excluding  yyx,  yyy 
and  yyz  components  which  never  exceed  0.1  x  10"’  esu. 
Tilt  greatly  modifies  the  pattern  of  components,  so  that 
for  example  even  at  20°  the  xzz  component  almost 
equals  the  zzz  component  in  magnitude,  with  the  zxx 
component  not  much  smaller.  Presumably  effects  of  this 
sort  correspond  to  the  noncentrosymmetry  in  the  plane 
indicated  by  the  SHG  measurements. 
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Abstract  Results  are  presented  on  the  design,  preparation  and  characterization 
of  Langmuir-Blodgett  films  formed  from  stearyl  (S)  derivatives  of  DAN,  MAP 
and  NPP.  Molecular  dynamics  modelling  of  trilayers  composed  of  simplified 
molecules  reveals  buckling  at  the  interface  owing  to  disordering  of  the  strongly 
interacting  head  groups.  Model  calculations  show  weak  dielectric  interactions 
between  neighbouring  layers,  and  reveal  major  effects  of  molecular  tilt  on  the 
quadratic  and  cubic  susceptibilities.  SDAN  forms  good  films  but  SNPP  and 
SMAP  do  not  form  stable  films  even  in  mixtures  with  stearic  or  behenic  acid. 
All  three  materials  give  SHG  as  powders  and  as  LB  films.  The  film  SHG 
coefficients  are  some  10%  of  those  for  crystals  of  the  underivatized  molecules. 


1.  INTRODUCTION 


Various  organic  molecules  such  as  4-iV,N-dimethylamino-3-acetamidonitrobenzene 
(DAN,  1),  2,4-dinitrophenyH-alanine  methyl  ester  (MAP,  2)  and  (5)-N-(4-nitro- 
phenyl)-prolinol  (NPP,  3)  are  known  to  give  crystals  with  high  nonlinear  optical 
response.  However,  crystal  growth  is  time-consuming  and  imposes  constraints  on  the 
possible  molecular  packings  and  lattice  symmetries.  Since  these  molecules  already 
possess  potentially  hydrophilic  groups,  we  have  adopted  the  strategy  of  derivatizing 
them  with  hydrophobic  stearyl  chains  to  produce  SDAN  (4),  SMAP  (5)  and  SNPP  (6). 


1:  DAN,  R  =  CH3 
4:  SDAN,  R  =  CnHj, 


NHCHMeCOR' 


NO2 

2:  MAP,  R'  =  CH, 

5:  SMAP,  R'  =  C„H370 


R"0 


3:  NPP,  R"  =  H 
6:  SNPP,  R"  =  C,7H3iCO 
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The  aim  is  to  deposit  the  derivatives  as  Langmuir-Blodgett  films  to  retain  the  high 
molecular  response  while  affording  more  control  of  the  material  structure. 

Our  approach  has  four  strands:  (i)  dynamic  modelling  of  films;  (ii)  calculation 
of  optical  properties;  (iii)  synthesis  of  film-forming  molecules  and  deposition  of  films; 
and  (iv)  linear  and  nonlinear  optical  characterization  of  films. 


2.  MOLECULAR  DYNAMICS  SIMULATIONS 

The  general  aim  of  the  simulations  is  to  understand  how  interactions  between 
molecules  control  the  structure  of  multilayer  LB  films.  Modelling  has  been  used  to 
characterize  the  two  types  of  interface,  between  chain  heads  and  between  the 
hydrophobic  chain  ends.  We  have  examined  structural  correlations  between  the  layers 
and  also  the  influence  of  structural  defects  of  the:  trilayer  structure. 

Our  approach  is  to  use  a  simplified  but  computationally  tractable  model  for  the 
amphiphiles  which  reflects  the  essential  molecular  features,  i.e.  connectivity,  van  der 
Waals  interactions,  and  a  realistic  representation  of  chain  flexibility.  In  this  initial 
study  the  chemically  complex  head  groups  are  represented  by  a  large  single  interaction 
site  with  parameters  chosen  to  reflect  the  bulkiness  and  intensity  of  the  associated 
interactions.  The  amphiphiles  consist  of  18  sites  linearly  connected  by  rigid  bonds 
(head  group  plus  16  methylene  groups  and  one  terminal  methyl  group).  Lennard-Jones 
12-6  potentials  are  used  throughout  with  a  cut-off  of  3.0  (in  units  of  the  appropriate 
cr).  The  parameters  for  the  head  group  interactions  are  =  500  K  is  the 
Boltzmann  constant),  =  4.611  nm,  and  =  0.100  kg  mol*',  and  for  the  tail  are 
ejk^  =  72  K  and  cr,  =  0.3923  nm.  The  usual  Berthelot  combining  rules  are  used  for 
cross  interactions.  Interactions  with  the  flat,  structureless  surface  are  modelled  using 
a  9  -  3  potential  with  the  same  parameters  as  used  previously:'  =  90  K  and  (Tbs 
=  0.36  nm,  except  that  the  head  groups  are  attracted  by  a  stronger  potential  ei^/k^  = 
302.132  K.  No  cut-off  is  used  for  the  surface  potential. 

The  method  of  preparation  is  to  build  and  equilibrate  the  film  layer  by  layer 
with  the  heads  of  the  first  layer  pointing  towards  the  surface.  Each  layer  contains  64 
chains  for  a  total  of  3456  particles.  The  chain  heads  are  initially  arranged  in  a 
hexagonal  lattice  with  the  rest  of  the  molecule  normal  to  the  surface  plane  and  in  the 
all-trans  state.  The  surface  area  per  head  group  is  23.2  AVchain.  Periodic  boundary 
conditions  are  imposed  in  the  X  and  Y  directions  and  the  minimum  image  convention 
is  employed.  The  first  16,000  time  steps  (40  ps)  are  discarded  for  equilibration.'  The 
data  are  obtained  from  100  ps  equilibrium  simulations. 

The  tilt  angle  d  is  determined  from  the  angle  between  the  surface  normal  and 
the  vector  between  the  fourth  and  fifteenth  beads.  At  300  K  its  mean  value  in  the  first 
layer  (closest  to  the  surface),  <5,  >  =  39°,  is  as  expected  from  the  simple  constant 
projected  area  model.^  The  smaller  values  of  <02>  =35°  and  <dj>  =  31°  in  the 
second  and  third  layers  indicate  that  each  layer  is  successively  more  disordered,  a  trend 
reflected  in  the  density  of  units  normal  to  the  surface  plane  G(Z)  presented  in  Figure 
1.  The  first  layer  is  highly  ordered  with  a  sharp  first  peak  corresponding  to  the  head 
groups.  Indeed,  seventeen  separate  peaks  corresponding  to  the  tail  groups  can  be 
resolved.  The  second  and  third  layers  are  less  ordered  although  the  head  peaks  are  still 
quite  distinct.  The  radial  distribution  function  projected  on  to  the  surface  plane  G(R) 
shows  that  the  hrads  in  the  first  layer  ^e  held  tightly  in  the  hexagonal  lattice  structure. 
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whereas  those  in  the -second  and  third  layers  show  almost  liquid-like  disorder. 
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FIGURE  1.  Density  profile  G(Z)  in  the  Z  direction  normal  to  the  LB  film. 

The  disordering  of  the  strongly  interacting  head  groups  can  be  related  to  the 
metastability  of  the  layer  structure  which  results  in  spontaneous  buckling  of  the 
interface.  TTiis  underlines  the  importance  of  incorporating  orientation-dependent  head 
group  int^ctions  in  order  to  stabilize  the  interfaces  better.  Some  indication  of  the 
extent  of  buckling  is  given  by  the  difference  in  half  widths  of  the  head  group  peaks  in 
(/(Z)  for  the  second  and  third  layers  as  compared  to  the  first  layer.  The  value  of  0.2  - 
0.3  nm  is  similar  in  magnitude  to  that  observed  experimentally  for  LB  films  by  atomic 
force  microscopy.’ 


3.  CALCULATION  OF  FILM  PROPERTIES 

Linear  and  nonlinear  optical  response  have  been  calculated  for  model  film  structures. 
'  The  molecules  are  represented  as  a  string  of  spherical  beads  to  model  their  size,  shape 
r  and  orientation,  the  number  of  beads  giving  the  axial  ratio  of  the  molecule.  The  long 
L  axes  are  initially  packed  in  a  triangular  array  normal  to  the  layers,  but  are  allowed  to 
i  tilt  by  an  angle  6  along  the  nearest-neighbour  direction,  as  suggested  by  the  molecular 
dynamics  calculations.  The  initial  structure  is  hexagonal  and  becomes  monoclinic  once 
C  the  molecules  tilt 

1  The  optical  properties  are  determined  by  the  molecular  properties,  transformed 
-  to  the  film  axis  system,  and  by  the  local  polarizing  electric  fields  in  the  layers.  The 
r  local  fields  depend  on  the  planewise  Loren  tz-factor  tensors  L(g),  which  are  sums  of 
?  dipole  interactions  between  a  molecule  in  layer  0  and  the  molecules  in  layer  g.* 
■r  Components  of  L(0)  (calculated  for  an_axial  lado  of  5)_refeQed  to  monoclinic  xyz 
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axes,  with  r-normal-to  the. layers  and  xz  the  plane  of  tilt,  are  shown  in  Table  I. 
Table  I.  Planewise  sum  components  L^(0)  for  various  tilts  6. 


e 

XX 

xz 

yy 

zz 

0° 

0.491 

0 

0.491 

-0.982 

20° 

0.440 

0.093 

0.467 

-0.907 

0 

0 

0.281 

0.148 

0.420 

-0.701 

At  zero  tilt  the  values  approach  (0.5,  0.5,  0),  which  is  the  depolarization  factor 
for  a  needle-shaped  Lorentz  cavity  appropriate  to  elongated  molecules.  (Results 
reported  previously^  contained  a  numerical  error.)  Tilt  about  y  affects  Lyy  little,  but 
induces  a  nonzero  and  brings  L„  and  closer.  For  g  =  1,  the  values  are  at  most 
1  %  of  the  corresponding  ones  in  Table  I,  and  so  can  be  neglected  in  practice. 

Optical  properties  are  calculated  using  a  model  molecular  polarizability  a  and 
first  and  second  hyperpolarizabilities  /3  and  7  having  the  same  axial  ratio  as  the 
molecules.  The  local  field  factors  are  1,0  to  1.2,  with  no  strong  dependence  on  tilt.  Of 
the  refractive  indices,  ity  varies  little,  while  increases  and  n.  decreases  with 
increasing  tilt,  reflecting  the  variation  of  the  polarizability  with  tilt.  Both  quadratic  and 
cubic  nonlinear  susceptibilities  have  been  calculated,  with  results  as  summarized  in 
Tables  II  and  m. 


Table  H.  Quadratic  susceptibility  components  Xc^pm  V*  for  various  tilts  B. 
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Table  m.  Cubic  susceptibility  components  Xc^(P^  V‘‘)^  for  various  tilts  6 
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The  pattern  of  components  varies  significantly  with  tilt. -This  reflects  predominantly 
the  effect  of  tilt  on  the  hyperpolarizability  components.  The  component  Xxxx  <ioes  not 
depend  monotonically  on  6  because  changes  sign  as  6  increases,  and  nor  does  Xzzzz 
because  the  direct  and  cascading  contributions  (which  depend  on  y  and  0^  respective¬ 
ly’)  have  opposite  dependences  on  6.  The  cascading  term  can  be  as  large  as  the  direct 
one,  and  hence  should  not  be  ignored  in  predicting  the  cubic  susceptibility. 


4.  SYNTHESIS  AND  FILM  DEPOSITION 

SDAN  (4)  is  synthesized  from  l-fluoro-5-nitroaniline  by  (i)  nucleophilic  displacement 
of  fluoride  in  ethanol  followed  by  (ii)  amidation  with  stearic  acid,  dicyclohexyl- 
carbodiimide  and  catalytic  4-pyrrolidinopyridine  in  dichloromethane.’  SMAP  (5)  is 
prepared  similarly  from  2,4-dinitro-fluorobenzene  by?  (i)  reaction  with  D-alanine  in 
ethanol  followed  by  (ii)  esterification  with  stearyl  alcohol.  SNPP  (6)  presents  more 
difficulty,  but  can  be  prepared  by  (i)  amination  of  4-fluoro-nitrobenzene  with  D-proline 
in  ethanol  followed  by  (ii)  esterification  with  stearic  acid.  All  three  materials  are 
crystalline  solids  at  room  temperature  and  give  satisfactory  NMR,  IR  and  micro- 
analytical  data  after  purification  by  column  chromatography  on  silica  gel. 

Dipping  experiments  are  carried  out  using  a  Joyce-Loebl  dual  A-B  trough. 
Well-ordered  homogeneous  non-centrosymmetric  multilayers  of  SDAN  are  deposited 
onto  silanized  glass  plates,  but  dipping  films  for  SMAP  and  SNPP  proves  more 
difficult  and  less  reliable.  The  films  are  typically  Y-type  for  SMAP  and  SN?"^  and  X- 
type  for  SDAN,  though  both  Y  and  Z-type  films  can  also  be  obtained  for  SDAN  by 
careful  control  of  the  experimental  conditions.  The  films  are  characterized  by  X-ray 
diffraction,  UV/visible  spectroscopy  and  optical  measurements,  which  together  clearly 
show  tilt  angles  ranging  from  near  zero  to  -30'’,  depending  on  the  dipping 
temperature  and  the  pH  of  the  subphase. 


Isotherm  of  SNPP 
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FIGURE  2.  Surface  isotherm  for  SNPP  indicating  monolayer  formation. 
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The  -highly  capricious  nature  of  the  films  formed  by  pure  SMAP  and  SNPP  led 
us  to  investigate  their  film-forming  properties  more  closely.  Systematic  variation  of 
temperature  and  pH  gives  some  isotherms  (e.g  at  8  ®C  and  pH  11:  Figure  2)  that 
indicate  the. formation 'qfca  fl  monolayer  for  SNPP,  but  no  conditions  could  be 

found  that  would  allow  routine  successful  dipping  of  films;  the  monolayer  seems  to 
collapse  over  a  period  of  about  20  minutes  at  reasonable  dipping  pressures.  Equimolar 
mixtures  of  SNPP  with  stearic  acid  or  behenic  acid  (which  sometimes  stabilize  films) 
deposit  better,  but  give  areas  per  molecule  suggestive  of  the  fatty  acid  alone,  a 
conclusion  supported  by  the  observation  of  SNPP  crystallites  apparently  squeezed  out 
of  the  film.  For  SMAP,  no  conditions  are  found  that  give  a  floating  monolayer. 
Equimolar  mixtures  with  behenic  acid  give  the  same  sort  of  behaviour  as  for  SNPP, 
indicating  expulsion  of  SMAP  from  the  floating  monolayer. 

The  differences  in  behaviour  of  these  three  p-nitroaniline  derivatives  suggest 
that  the  position  of  attachment  of  the  long  alkyl  chaineto  the  head  group  controls  the 
conformation  and  arrangement  of  the  molecules  in  a  manner  that  is  critical  for  LB  film 
formation.  We  are  currently  investigating  how  to  improve  the  stability  and  reproduc¬ 
ibility  of  the  films  by  methods  including  (i)  varying  the  position  of  attachment  of  the 
tails  to  the  NLO-active  head  groups,  (ii)  using  multiple  tails,  and  (iii)  organising  the 
head  groups  through  specific  directional  and  recognitive  interactions. 


5.  FILM  CHARACTERIZATION 


Second-harmonic  generation  (Figure  3)  is  a  sensitive  probe  of  structure  within  the  LB 
nonlinear  optical  films.  A  Nd:YAG  laser  is  used  which  emits  Q-switched  pulses  of  1 
mj  energy  in  10  ns  at  a  wavelength  of  1064  nm  with  a  repetition  rate  of  30  Hz.  The 
pulses  have  a  Gaussian  spatial  profile,  with  stability  of  the  order  of  1%.  The  light  is 
p-polaiized  and  focused  onto  the  sample  with  a  spot  size  of  300  micron.  The  angle  of 
incidence  onto  the  sample  is  stepper-motor  controlled.  After  passage  through  an 
infrared-blocking  KG3  low  pass  filter,  neutral  density  filters  as  required  and  a  1  nm 
bandwidth  interference  filter  for  the  second  harmonic  at  532  nm,  the  harmonic  signal 
is  detected  by  a  photomultiplier.  A  boxcar  integrator  is  used  to  gate  the  detection 
system  and  reduce  noise. 
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FIGURE  3.  Experimental  arrangement  for  studying  SHG. 
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The  powders  of  SDAN,  SMAP  and  SNPP  still  exhibit  SHG,  and  so  the  stearyl 
group  does  not  promote  centrosymmetry,  in  contrast  to  some  results  for  other  p-nitro- 
aniline  derivatives.  The  LB  films  of  SDAN,  SMAP  and  SNPP  also  exhibit  SHG. 
However,  since  SDAN  produces  the  best  films,  it  has  been  the  subject  of  the  major 
part  of  this  study.  The  dependence  of  SHG  on  the  angle  of  incidence  onto  an  LB  film 
can  be  analysed  to  find  the  orientation^  the  nonlinear  chromophores  within  the  layers. 
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SHG  results transmission  through  an  SDAN  film  coated  on  both  sides  of  the 
substrate  are  shown  in  Figure  4.  The  interference  between  the  second  harmonic 
produced  at  the  ^o  surfaces  creates  fringes  with  varying  angle  of  incidence.  The 
significant  harmonic^sagnal  at  normal  incidence  indicates  non-centrosymmetric 
molecular  packing  in  the  plane  of  the  film,  and  the  angular  dependence  of  the  SHG 
efficiency  indicates  that  the  nonlinear  chromophores  lie  precisely  in  the  plane  of  the 
film  at  the  interface  of  each  layer.  The  nonlinear  axis  of  the  SDAN  molecule  is  not 
parallel  to  the  molecular  axis  defined  by  the  hydrocarbon  chain.  The  measured 
intensity  does  not  vary  across  the  film  perpendicular  to  the  dipping  direction,  or  along 
the  film  parallel  to  it,  to  within  the  accuracy  of  the  measurements. 
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FIGURE  4.  Interference  fringes  from  SHG  in  an  SDAN  LB  film. 

LB  films  of  SDAN  deposited  on  a  single  trough  exhibit  SHG  at  normal 
incidence  similar  to  the  materi^  DCANP,*  also  a  derivative  of  p-nitroaniline.  The 
:  nonlinear  coefficient  of  SDAN  is  d  =  7  x  10*’  esu,  very  similar  to  the  value  of 
8  X  10*’  esu  reported  for  DCANP.  For  SDAN  the  SHG  coefficient  is  about  10%  of 
I  that  of  the  DAN  crystal,  without  the  stearyl  chain.  For  SNPP  the  coefficient  relative 
I-  to  NPP  is  slightly  greater  than  10%.  Preliminary  measurements  on  SMAP  have  shown 
r  high  SHG  conversion  efficiencies,  but  with  evidence  of  microcrystallite  formation. 

The  films  of  SDAN  have  a  very  high  optical  damage  threshold,  in  excess  of  1 
'  GW  cm*^.  They  have  excellent  environmental  stability,  with  no  detectable  deterioration 
-of  nonlinear  response  after  over  a  year  of  storage  at  room  temperature  including 
C  exposure  to  fluorescent  lighting.  Systematic  study  of  the  film  quality  as  a  function  of 
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deposition  conditions  has  not  yet  allowed  the  construction  of  multilayer  stacks  over  50 
layers  thick.  The  films  have  a  slightly  milky  appearance  due  to  the  formation  of 
Rayleigh  scatter  domains  within  the  film  structure.  Waveguiding  experiments  have 
therefore  not  yet  been  possible. 


6.  CONCLUSIONS 

Molecular  dynamics  simulations  provide  insight  into  the  structure  and  perfection  of  LB 
films  and  how  these  relate  to  the  molecular  structure.  Model  calculations  show  how 
the  molecular  nonlinear  response  expresses  itself  in  the  films,  with  particular  reference 
to  the  effect  of  tilt.  Molecules  that  produce  SHG-active  crystals  have  been  deiivatized 
to  produce  LB  films,  though  these  are  not  always  stable.  The  films  have  SHG 
coefficients  that  are  10%  of  those  for  the  crystals,  with  good  damage  resistance  and 
environmental  stability. 

Current  work  is  aimed  at  studying  in  more  detail  how  molecular  structure 
affects  film  NLO  response,  in  order  to  improve  molecular  design,  and  at  designing  and 
synthesizing  new  species  to  give  uniform  thick  multilayer  films  with  better  optical 
performance. 
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The  sample  size  dependence  and  time  convergence  of  property  values  of  simple  hydrocarbon 
monolayers  has  been  investigated  by  performing  molecular  dynamics  simulations  at  four 
surface  coverages  for  systems  ranging  in  size  from  16  to  1(X)  chains  and  for  trajectories  as  long 
as  200  ps.  Detailed  studies  of  the  tilt  angle  indicate  that  systems  with  64  chains  simulated  for 
40  ps  (in  addition  to  a  40  ps  equilibration  stage)  are  large  enough  to  obtain  statistically 
meaningful  results.  The  equilibrium  tilt  angle  is  slightly  sensitive  to  the  boundary  conditions. 


The  investigation  of  the  structure  of  Langmuir  and 
Langmuir-Blodgett  films  has  now  become  an  extremely  ac¬ 
tive  area  of  research'  driven  by  potential  applications  in  the 
fabrication  of  electronic  and  electro-optic  devices.  Langmuir 
films  are  molecular  monolayers  at  the  liquid-gas  interface 
(usually  water-air)  whereas  Langmuir-Blodgett  films  are 
mono  or  multilayer  films  deposited  from  the  liquid-air  inter¬ 
face  onto  a  solid  substrate  such  as  silver-plated  mica.  In  or¬ 
der  to  form  a  Langmuir  film  a  molecule  must  be  amphiphilic 
so  that  its  hydrophilic  “head”  is  immersed  into  the  water 
phase  while  its  hydrophobic  tail  points  into  the  gas  phase. 
The  classic  materials  used  to  make  Langmuir  films  are  the 
fatty  acids  such  as  arachidic  acid,  CH3(CH2)|«C02H,  and 
its  metal  salts. 

Several  simulations  of  idealized  monolayers  have  been 
reported  in  recent  years. Each  research  group  has,  how¬ 
ever,  used  a  different  molecular  model,  making  it  difficult  to 
compare  results  and  to  judge  the  correctness  of  the  conclu¬ 
sions.  In  order  to  assess  the  influence  of  sample  size,  initial 
lattice  type  and  the  convergence  of  property  values  with  time 
we  have  carried  out  a  systematic  examination  of  one  of  these 
published  models^  *  *  for  N=20  unit  hydrocarbon  mono- 
layers. 

In  this  “united  atom”  model  the  hydrogens  are  col¬ 
lapsed  into  the  carbons  so  that,  for  example,  a  methyl  group 
is  replaced  by  a  spherical  bead  of  mass  0.01 5  kg  mol  ~  '  and  a 
methylene  group  by  a  spherical  bead  of  mass  0.014 
kg  mol  ~  '.  The  connectivity  of  the  chains  is  maintained  by 
joining  the  interaction  sites  by  N-l  rigid  bonds  (kept  at  a 
constant  magnitude  of  0. 153  nm,  with  a  relative  tolerance  of 
10  “  by  the  SHAKE*’  procedure  for  constraint  dynam¬ 
ics). 

Each  chain  has  N-2  bond  angles,  0,  ,  which  are  subject  to 
the  following  harmonic  potential 

U(di)=o.5kg{e,-0„)\  (1) 

Here,  kg  and  are  the  equilibrium  bond  angle  force  con¬ 
stant  and  bond  angle,  respectively.  These  have  the  values'" 
kg  =  520  kJ  mol  "  '  rad  ~  ^  and  =  109.4T.  A  torsional  po- 
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tential  (7(cr )  for  the  dihedral  angle  a  between  the  two  planes 
formed  by  four  consecutive  chain  units  is  included. ' ' 

f/(a)  =  Co  -f-  c,  cos(a)  -t-  c,  cos-(a) 

-I-  Cj  cos^(a)  -b  Cg  cos’(£r)  -(-  Cj  cos*(a),  (2) 

where  Cq,  c„  C2,  Cj,  Cg,  and  c,  have  the  values  9.2789, 12. 1557, 
-  13.1201,  -  3.0597,  26.2403,  and  -  31.4950  in  units  of 
kJ  mol  ■ ',  respectively. 

In  addition  to  the  above  potential  interactions  all  pairs 
of  units  on  a  given  chain  which  are  separated  by  at  least  four 
bonds  and  all  pairs  on  different  chains  interact  through  a 
Lennard-Jones  12-6  potential  which  is  truncated"  at  2.294 
O’.  The  values  for  the  Lennard-Jones  constants"  are 
€/kg  =  72  K  ( is  Boltzmanns’  constant)  and  o  =  0.3923 
nm. 

In  the  modeP'"’"  all  of  the  units  interact  with  a  flat,  struc¬ 
tureless  surface  (set  at  Z  =  0)  by  a  functional  form  of  the 
Lennard-Jones  potential  which  results  from  integration  over 
the  surface'^  (viz.  9-3  instead  of  the  usual  12-6). 

Ugg  =  20€gs[(,(Tgg/Z)'* — (o’ns/Z)^].  (3) 

We  have  selected  the  same  set  of  parameter  values  as  em¬ 
ployed  by  Bareman  and  Klein"  Css^kg  =  90  K  and 
ags  =  0.36  nm.  No  cutoff  is  used  for  the  surface  potential. 

The  systems  initially  have  their  chain  “heads”  arranged 
in  either  a  triangular  or  a  square  lattice  with  the  rest  of  the 
molecule  normal  to  the  Z  =  0  plane  and  in  the  all  trans  state. 
The  Verlet  algorithm'^  with  a  time  step  of  2.5  fs  was  used  to 
integrate  Newton’s  equations  of  motion.  Periodic  boundary 
conditions  were  imposed  in  the  X  and  Y  directions  and  the 
minimum  image  convention  was  employed.  Following  Bare¬ 
man  et  al.*  the  first  16  (XX)  time  steps  (40  ps)  have  been 
discarded  for  equilibration.  Molecular  dynamics  simula¬ 
tions  for  80  ps  on  a  single  processor  of  a  Silicon  Graphics  260 
GTX  require  about  four  days  of  computer  time  for  64  chain 
systems. 

The  surface  coverage  is  conveniently  defined  by  the  area 
per  chain  .  We  have  carried  out  simulations  atA^  values  oP 
19.6, 20.8, 22.0,  and  23.2  A’/chain  for  sample  sizes  rang¬ 
ing  from  1 6  to  ICX)  chains.  For  the  two  largest  values  of  xl,  we 
have  used  200  ps  simulations  which  appear  to  be  ten  times 
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FIG.  I.  The  time  variation  of  the  tilt  angle  for  a  64  chain  system  with 
=  22.0  A^/chain.  The  horizontal  line  is  the  final  average. 

longer  than  those  reported  by  Bareman  et  ai*  The  main 
quantities  which  have  been  monitored  are  the  temperature 
(T), the variousenergies  ( bondangle (EB ) , torsional  (ET ) , 
Lennard- Jones  {JsU  ),  surface  (ES ),  and  total  potential  en¬ 
ergy  (EP ) ),  the  percentage  of  trans  states  (conformations 
with  |o|<6(J’)  {%T),  and  the  tilt  angle  from  the  normal 
direction  {0).  The  tilt  angle  has  been  determined  from  the 


TABLE  I.  The  simulation  results  starting  with  a  triangular  lattice. 


=  19.6  (20  ps) 

16 

64 

(D/K 

<£I-/)/kJitiol-' 

(ESI/kJmol' 

<%r> 

<^)/deg 

302.31  ±0.95 
-  158.80  ±0.11 
-  17.09  ±0.33 

99.43  ±  0.32 

20.11  ±0.12 

300.30  ±0.45 
-  159.17  ±0.04 
-  16.70  ±  0.02 

99.39  ±0.02 

20.34  ±  0.04 

=  20.8  (  20  ps) 

16 

64 

<r>/K 

(EU)/\axniA' 

<£S)/kJmo|-' 

(%T) 

<^)/deg 

301.63  ±0.97 
-  154.88  ±0.11 
-  17.79  ±0.05 

98.05  ±  0.08 

26.06  ±  0.09 

297.75  ±  0.50 
-  156.66  ±0.06 
-  17.59  ±  0.02 

98.82  ±  0.02 

26.99  ±  0.05 

.4,  »22.0(200ps) 
16  64 

100 

(r)/ic 

(ES)AJmtA-' 

<%n 

<B>/aeg 

306.66  ±  1.05  310.09  ±0.53 

-  152.69  ±0.12  -  155.06  ±0.05  - 

-18.58  ±0.04  -18.91  ±0.02 

98.32  ±  0.06  98.44  ±  0.03 

31.05  ±  0.08  32.85  ±0.04 

306.25  ±0.36 

-  155.05  ±0.04 

-  18.86  ±0.02 
98.49  ±0.02 
3178  ±0.03 

23.2  (200 ps) 
16  64 

100 

<T)/K 

<£L/)/iLf  mol  - ' 
<£5)/U  mol  ' 

<%r) 

<4>/«lei 

302.61  ±  1.06  300.39  ±0.54 

-  147.47  ±  0.12  -  153.29  ±0.06  - 
-  19.33  ±  0.05  -  19.90  ±  0.02 

97.12  ±0.09  98.28  ±004 

33.96  ±  0.08  36.89  ±  0.04 

305.83  ±0.41 
-  153.85  ±0.04 
-20.01  ±0.02 
98.36  ±0.03 
37.31  ±0.04 

J.  Ch«m.  Phys.,  Vof. 
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FIG.  2.  The  extrapolation  of  the  tilt  angle  for  different  areas  per  chain  for 
the  triangular  lattice. 


angle  the  vector  between  the  fourth  and  seventeenth  beads 
makes  with  the  normal  direction.*  The  chain  ends  are  not 
used  in  the  calculation  of  the  tilt  angle  since  the  probability 
for  gauche  configurations  at  the  ends  can  be  quite  high.^"’  ’* 
The  angular  brackets  indicate  that  the  value  of  the  system 
property  is  time  averaged  over  the  equilibrated  molecular 
dynamics  trajectory.  The  variation  of  tilt  angle  with  time  is 
shown  in  Fig.  1  for  the  2(X)  ps  following  the  initial  40  ps 
equilibration  phase  in  a  64  chain  system  with  =  22.0 
AVchain.  This  figure  suggests  that  the  value  of  32.85  indi¬ 
cated  by  the  horizontal  line  is  close  to  the  equilibrium  aver- 
a^. 

The  simulation  results  for  the  triangular  lattice  are  pre¬ 
sented  in  Table  I.  The  errors  indicated  are  one  standard  devi¬ 
ation  of  the  mean.  Data  were  gathered  at  intervals  which  are 
much  larger  than  the  correlation  time  so  that  it  is  sufficient 
to  use  statistical  methods  for  independent  quantities.  The 
error  bars  on  the  larger  chain  system  values  are  smaller  than 
those  for  the  smaller  chain  values  as  is  expected  from  statisti¬ 
cal  considerations.  For  all  of  the  results  with  16  chains  the 
range  of  the  potential  in  the  smaller  of  the  cell  dimensions  is 
slightly  less  than  the  Lennard-Jones  cutoff;  we  do  not  expect 
this  to  have  a  strong  influence  on  the  observed  properties. 

The  results  show  that  16,  64,  and  100  chain  systems 
display  approximately  the  same  energy  values  per  mole.The 
internal  chain  energies  such  as  the  bond  angle  energy  and  the 
torsional  energy  are  hardly  affected  by  the  change  in  but 
the  Lennard-Jones  energy  becomes  less  negative  and  the  sur¬ 
face  energy  becomes  more  negative  as  A^  increases.  These 


TABLE  II.  The  tilt  angle  at  dilTeient  areas  per  chain. 


Extrapolation 

Eq.  (4) 

Bareman  and  Klein 
(Ref.  6) 

19.6  AVchain 

20.42  ±0.07 

19.98 

=20 

20.8  AVchain 

27.30  ±0.07 

27.68 

=27  * 

22.0  AVchain 

33.20  ±0.03 

33.15 

=  31 

23.2  A^/chain 

37.92  ±  0.04 

37.44 
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0  4  8  12  16  20  0  40  00  120  160  200 

TIME  /  PS  TinE  /  PS 

FIG.  3.  Time  variation  of  the  running  average  of  the  tilt  angle  for  different 
areas  per  chain;  solid  line — triangular  lattice,  broken  line — square  lattice. 
The  horizontal  lines  are  the  final  average  values. 


observations  reflect  the  enhanced  influence  of  the  surface  as 
increases.  Similarly ,  the  average  tilt  angle  increases  as 
increases  because  the  chains  have  more  room  in  which  to 
cant  over.  The  percentage  of  gauche  states  increases  slightly 
from  about  1%  to  1.5%  as  increases  since  more  space  is 
now  available  for  gauche  defects  to  develop.  The  percentage 
of  trans  states  found  in  these  molecular  dynamics  simula¬ 
tions  is  in  good  agreement  with  other  workers.  ^  ‘ 

Since  the  tilt  angle  is  a  crucial  property  of  film  systems 
we  have  investigated  its  behavior  in  more  detail.  The  tilt 
angle  values  for  finite  have  been  used  to  obtain  the  tilt 
angle  values  for  infinite  by  fitting  a  weighted  least- 
squares'^  linear  line  in  l/iV,  to  each  set  of  data  in  Table  I. 
The  resulting  intercepts  are  20.42  ±  0.07,  27.30  ±  0.07, 
33.20  ±  0.03,  and  37.92  ±  0.04  for  A,  =  19.6, 20.8, 22.0  and 
23.2  A^/chain,  respectively.  This  extrapolation  procedure  is 
illustrated  in  Fig.  2.  The  fits  at  the  higher  area  per  chain 
values  indicate  that  for  a  small  number  of  chains  (.N^  <7) 
the  tilt  angle  of  the  23.2  A^/chain  system  will  be  smaller  than 
that  of  22.0  A  Vchain  system.  This  reflects  the  large  influ¬ 
ence  of  the  boundary  conditions  in  small  chain  systems. 

Bareman  and  Klein''  suggest  that  the  variation  of  the  tilt 
angle  with  is  given  by 

{0}=cos-'(A^/A,),  (4) 

where  A^  is  the  area  per  chain  at  maximum  thickness  (i.e., 
when  all  the  molecules  are  normal  to  the  surface).  They 


found  for  the  “united  atom”  model  that  A^  =  18.42 
A  Vchain.  Table  II  lists  the  extrapolated  molecular  dynam¬ 
ics  data,  the  predictions  of  Eq.  (4)  and  the  molecular  dy¬ 
namics  results  of  Bareman  and  Klein*  for  =  90;  their 
values  have  been  read  from  their  Fig.  3.  The  values  predicted 
by  Eq.  (4)  are  in  excellent  agreement  with  our  extrapolated 
molecular  dynamics  data.  The  value  at  22.0  A^/chain  re¬ 
ported  by  Bareman  and  Klein*  ( ~  3 1 )  is  significantly  lower 
than  that  found  in  this  work  or  predicted  by  Eq.  (4).  We 
suspect  that  their  trajectory  was  not  long  enough  adequately 
to  sample  the  tilt  angle  at  this  area  per  chain. 

Since  the  tilt  angle  results  for  the  64  chains  are  in  excel¬ 
lent  agreement  with  Bareman  and  Klein’s*  calculation  using 
90  chains  ( at  the  two  lowest  areas  per  chain )  and  are  also  in 
excellent  agreement  with  the  100  chain  systems  (at  the  two 
largest  areas  per  chain)  we  can  conclude  that  64  chain  sys¬ 
tems  are  large  enough  to  give  accurate  calculations  of  the  tilt 
angle. 

The  convergence  of  the  tilt  angle  with  time  is  illustrated 
in  Fig.  3  for  the  64  chain  systems.  Even  at  the  highest  area 
per  chain  of  23.2  AVchain  the  running  average  has  con¬ 
verged  by  40  ps.  Thus,  it  is  sufficient  to  simulate  64  chain 
systems  for  40  ps  ( in  addition  to  a  40  ps  equilibration  stage ) 
in  order  to  have  good  statistics  on  the  tilt  angle. 

The  simulations  have  been  repeated  for  64  chain  sys¬ 
tems  starting  with  the  heads  initially  arranged  in  a  periodic 
square  lattice  and  the  results  are  contained  in  Table  III.  At 
the  two  lower  densities  the  runs  were  performed  for  40  ps  ( in 
addition  to  a  40  ps  eqfilibration  stage).  At  the  higher  densi¬ 
ties  results  from  triangular  lattice  simulations  suggest  that 
20  ps  is  sufficient  to  accumulate  satisfactory  statistics.  The 
convergence  of  tilt  angle  with  time  is  included  in  Fig.  3.  It  is 
seen  that  convergence  is  achieved  for  all  the  areas  per  chain 
studied.  When  one  compares  Tables  I  and  III  one  notes  that 
(£S )  and  {p )  for  the  square  lattice  is  always  larger  than  for 
the  triangular  lattice  at  the  same  values  of  A  ^  (the  difference 
in  (^  >  increases  from  about  0.4”  to  1 .2’  as  decreases  from 
23.2  to  19.6  AVchain). 

That  is  at  first  somewhat  surprising  since  the  triangular 
lattice  is  the  closed-packed  arrangement  which  for  a  given 
coverage  might  be  expected  to  allow  more  efficient  packing 
of  the  tails  and  hence  lead  to  greater  tilt.  Although  the  mole¬ 
cules  are  free  to  move  in  the  surface  plane,  for  the  square 
lattice  simulations  they  are  prevented  from  taking  up  a  close- 
packed  arrangement  by  the  boundary  conditions.  Neverthe¬ 
less  enhanced  tilt  can  still  be  achieved  along  the  diagonals  of 
a  square  lattice. 

The  research  reported  herein  has  been  sponsored  in  part 


TABLE  III.  The  simulation  results  for  a  64  chain  system  starting  with  a  square  lattice. 


19.6  (20  ps) 

20.8  (20  ps) 

22.0  (40  ps) 

23.2  (40  ps) 

<r>/K 

303.90  ±  0.35 

302.95  ±  0.49 

307.83  ±0.33 

302. 10  ±0.54 

(EU)AOmtH  ' 

-  159.31  ±0.03 

-  157.86  ±0.06 

-  155.77  ±0.05 

-  153.21  ±0.06 

(£S>/Idmol  ' 

-  16.91  ±0.02 

-  18.00  ±0.02 

-  19.01  -1-0.02 

-  19.93  ±0.02 

<%T) 

99.17  ±0.02 

98  69  ±  0.02 

98.82  ±0.03 

98.27  ±  0.03 

21.50  ±0.04 

28.30  ±0.05 

33.31  ±0.05 

37  29  ±0.04 
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CALCULATIONS  OF  NONLINEAR  OPTICAL 
PROPERTIES  OF  MODEL  LANGMUIR-BLODGETT  FILMS 


R.W.  Munn  and  M.M.  Shabat 

Department  of  Chemistry  and  Centre  for  Electronic  Materials, 

UMIST,  Manchester  M60  IQD,  U,K. 

ABSTRACT 

Calculations  have  been  performed  on  models  of  Langmuir-Blodgett  LB  films 
to  explore  how  molecular  nonlinear  response  expresses  itself  in  film  response.  The 
molecules  are  treated  as  a  string  of  s  beads  initially  arranged  normal  to  the  film  with 
their  long  axes  forming  a  close-packed  hexagonal  array,  but  they  are  then  allowed  to 
tilt  away  from  the  vertical.  The  molecular  polarizability  a  and  first  hyperpolarizability 
P  are  represented  by  one-dimensional  models  in  which  the  axial  components  are  s 
times  the  perpendicular  ones.  Because  of  the  layer  structure  of  the  films,  the  dipolar 
interactions  are  calculated  as  planewise  sums  between  one  molecule  and  all  the 
molecules  in  a  given  layer.  The  sums  have  been  calculated  for  various  values  of  s  as 
a  function  of  tilt;  they  are  negligible  except  for  interactions  of  a  molecule  with  its 
own  layer.  They  have  then  been  used  to  calculate  the  local  electric  field,  refractive 
index  and  quadratic  susceptibility  for  thin  films.  The  results  show  that  tilt  has  a 
major  effect  on  the  pattern  of  components  of  and  hence  on  measured  nonlinear 
optical  behaviour.  Such  information  can  be  used  to  help  guide  the  design  of  molecules 
to  achieve  desired  film  properties. 

1.  INTRODUCTION 

Langmuir-Blodgett  films  promise  a  number  of  advantages  for  applications  in 
molecular  electronics.  They  are  compact  and  highly  ordered,  and  their  thickness  and 
structure  can  be  controlled  by  depositing  the  desired  number  and  sequence  of 
monomolecular  layers.  However,  realizing  their  advantages  presents  the  challenge  of 
designing  and  synthesizing  molecules  which  impart  the  requisite  optical  or  electrical 
activity  while  affording  stable  high-quality  films. 

We  are  undertaking  a  coordinated  project  to  develop  LB  films  for  nonlinear 
optics.'  The  overall  strategy  is  to  take  molecules  which  give  crystals  with  high  NLO 
activity  and  to  derivatize  them  so  that  they  give  good  LB  films  which  retain  high 
NLO  activity.  Our  approach  has  four  strands; 

•  Molecular  dynamics  simulations  to  relate  film  structure  and  stability  to 
molecular  structure. 

•  Calculation  of  film  properties  from  film  structure  to  predict  NLO  activity. 

•  Synthesis  and  deposition  of  LB  films  to  evaluate  molecules  and  techniques. 

•  NLO  characterization  of  LB  films  to  assess  film  performance  and  device 
potential. 
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Eventually  the  aim  is  that  the  strands  should  operate  in  this  sequence,  with  each 
providing  feedback  to  the  earlier  ones,  but  this  is  not  yet  fully  achieved  while  the 
individual  strands  are  being  developed. 

The  present  paper  describes  work  on  the  second  strand.  \N'e  have  performed 
calculations  on  simple  models  of  LB  films  to  explore  how  molecular  response 
expresses  itself  in  film  response.  This  depends  on  two  factors:  the  molecular 
arrangements  and  the  molecular  interactions.^  We  therefore  first  describe  our  model 
for  film  structure.  Then  we  calculate  the  dipolar  interactions  within  the  film.  This 
finally  allows  us  to  calculate  the  film  properties,  with  particular  reference  to  the  effect 
of  molecular  tilt  away  from  the  normal  to  the  film. 

2.  MODEL 

Our  model  is  designed  to  give  the  film  structure  directly  from  the  molecular 
structure.  Molecules  which  yield  LB  films  are  typically  elongated,  and  as  such 
resemble  those  which  yield  liquid  crystals.  The  effect  of  molecular  elongation  has 
been  examined  previously  in  the  context  of  linear  optics  of  liquid  crystals^  and  later 
in  nonlinear  optics  of  molecular  crystals."*  The  present  model  derives  from  this 
earlier  work. 


Elevation 

Figure  1.  Sketch  of  the  model  LB  film  structure  with  molecules  vertical  for  i  =  5. 

Each  molecule  is  treated  as  a  string  of  s  spherical  "beads",  or  submolecules. 
In  the  basic  structure  the  molecules  in  each  layer  form  a  triangular  close-packed 
arrangement  with  their  long  axes  normal  to  the  layer.  Alternate  layers  are  centred 
over  the  triangle  formed  by  the  previous  layer  in  such  a  way  that  the  overall  structure 
has  hexagonal  symmetry.  For  simplicity,  each  layer  lies  entirely  above  the  preceding 
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one  rather  than  packing  into  the  recesses  of  the  triangles  in  the  layer  below,  so  that 
the  structure  does  not  reduce  to  hexagonal  close  packing  for  i  =  1 .  This  structure 
corresponds  to  a  crystal  with  lattice  parameters  a  =  b  =  2D,  c  =  2sD,  a  =  /3  = 
90°  and  y=  120°,  where  D  is  the  diameter  of  a  bead.  There  are  two  molecules  k  = 
1,  2  in  the  unit  cell,  with  submolecules  j  located  at  fractional  coordinates  [0,  0, 
(j  -  l)/2s])  for  ^  =  1  and  [1/3,  2/3,  1/2  +  (j  -  l)/2s)]  for  k  =  2,  where  j  = 

This  structure  is  illustrated  in  Figure  1  on  the  previous  page. 

Experiments  indicate,  however,  that  molecules  are  typically  tilted  away  from 
the  vertical,  and  molecular  modelling  confirms  this  tendency,  explaining  it  as  arising 
from  the  attraction  of  the  elongated  molecules  for  the  substrate.  The  resultant 
structure  depends  on  the  plane  in  which  this  tilt  occurs.  All  distinct  tilted  structures 
are  encompassed  by  tilts  in  planes  lying  between  the  nearest-neighbour  direction  (AB 
in  Figure  1)  and  the  next-nearest-neighbour  direction  (AC  in  Figure  1).  The  molecular 
modelling  indicates  that  next-nearest-neighbour  tilt  is  preferred,  and  we  restrict 
ourselves  to  such  structures. 

Now  LB  films  are  classified  as  X,  Y  or  Z-type.*  X  and  Z-type  films  deposit 
on  the  substrate  only  on  the  upstroke  or  downstroke  of  the  substrate  through  the 
floating  Langmuir  film,  while  the  commonest  (and  most  stable)  Y-type  films  deposit 
on  both  strokes.  It  is  often  assumed  that  Y-type  films  are  centrosymmetric  for  even 
numbers  of  layers  and  hence  cannot  exhibit  quadratic  nonlinearities.  This  is  correct 
for  directions  normal  to  the  film,  but  not  for  directions  parallel  to  the  film  when  the 
molecules  are  tilted.*  This  applies  whether  the  molecules  are  all  tilted  in  the  same 
sense,  to  alternate  in  the  direction  of  tilt,  in  which  case  the  orthogonal  symmetry  is 
retained.  (See  Figure  2).  We  consider  here  only  the  case  of  equal  tilt  0  away  from  the 
vertical.  The  lattice  parameters  are  then  a  -  b  =  2D,  c  =  2D  cos  6,  a  = 
cos  '[(\/3/2)  sin  S],  =  90°,  7  =  120°,  with  fractional  coordinates  unchanged. 


Tilt  angles  the  same 


Tilt  angles  alternating 


Figure  2.  Model  LB  film  structures  with  molecules  tilted  for  j  =  5. 
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We  also  need  a  model  for  the  molecular  response.  Calculations  of  molecular 
polarizability  and  hyperpolarizability  as  a  function  of  elongation  mostly  refer  to 
conjugated  molecules  as  opposed  to  the  molecules  envisaged  here  with  a  conjugated 
NLO  head  group  and  a  saturated  hydrocarbon  tail.  For  conjugated  molecules  the  axial 
molecular  response  is  typically  found  to  vary  superlinearly  with  the  molecular  length. 
For  illustrative  purposes  we  take  the  axial  response  in  our  molecules  to  be  simply 
proportional  to  the  length.  In  molecular  axes  X,  Y,  Z  the  polarizability  tensor  a  is 
then  diagonal,  with  components  (a,  a,  sa),  and  the  hyperpolarizability  tensor  0  has 
non-zero  components  $2ab  diagonal  in  AB,  with  components  (/?,  -  Sj3).  Though 

clearly  oversimplified,  these  assumptions  yield  an  axial  molecular  response  per  unit 
volume  which  is  independent  of  molecular  length,  so  that  changes  in  susceptibility  are 
directly  attributable  to  changes  in  structure. 

3.  PLANEWISE  DIPOLE  SUMS 

Linear  and  nonlinear  response  of  molecular  materials  depend  on  the  dipolar 
interactions  in  the  structure.’  These  relate  the  local  polarizing  field  at  a  molecule  to 
the  macroscopic  electric  field  E  which  appears  in  Maxwell’s  equations.  In  a  layered 
structure  it  is  natural  to  express  these  interactions  as  a  sum  of  planewise- dipolar 
interactions  giving  the  field  at  a  molecule  in  one  layer  due  to  dipoles  on  the  molecules 
in  the  same  or  another  layer,  calculations  for  layers  of  point  molecules  in  aromatic 
hydrocarbon  crystals  show  that  such  sums  fall  off  exponentially  with  distance  between 
the  molecule  and  the  layer  of  dipoles.*  This  gives  rise  to  the  idea  of  a  range  r  beyond 
which  the  interactions  are  negligible,  so  that  r  defines  the  number  of  coupled  planes 
to  be  included  in  treating  the  material  properties.’ 

We  have  calculated  planewise  sums  using  standard  expressions. Sums  are 
calculated  between  all  pairs  of  submolecules  and  then  averaged  to  give  sums  which 
reflect  the  molecular  size,  shape  and  orientation.  For  illustrative  purposes  we  report 
the  results  for  r  =  5  and  explore  how  the  sums  depend  on  the  tilt  6.  The  results, 
made  dimensionless  by  a  factor  v/4ir,  with  v  the  molecular  volume,  are  given  in 
Table  I. 


8 

a/? 

XX 

xz 

yy 

zz 

0“ 

0.491 

0 

0.491 

-0.982 

n 

o 

0 

0.429 

-0.099 

0.456 

-0.886 

o 

O 

0.284 

-0.184 

0.412 

-0.696 

Table  L  Planewise  dipole  tensor  sum  components  T^(g),  averaged  over  submolecules 
as  a  function  of  tilt  6  for  interactions  in  the  same  layer  (g  =  0),  expressed  relative 
to  cartesian  x  y  z  axes,  with  z  the  normal  to  the  layers  and  xz  the  plane  of  tilt. 
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The  diagonal  components  of  the  planewise  dipole  sums  satisfy  Tr  T(^)  =  0. 
They  are  markedly  anisotropic,  but  tilt  reduces  the  anisotropy  and  induces  nonzero 
omponents  T^(£),  consistent  with  monoclinic  symmetry.  For  elongated  molecules 
.ike  these,  the  appropriate  Lorentz  cavity  would  be  needle  shaped,  and  then  we  would 
have  the  xt  components  equal  to  1/2  and  the  zz  component  equal  to  -1;  this  is  clearly 
a  reasonable  approximation  to  the  results  for  zero  tilt.  We  find  that  the  components 
for  g  =  1  are  much  smaller  than  those  for  g  =  0  and  so  can  be  neglected  (i.e.  r  = 
1).  Results  for  values  of  s  up  to  10  are  very  similar.  We  therefore  conclude  that  the 
dielectric  properties  of  LB  films  are  dominated  by  interactions  within  one  layer,  with 
negligible  contributions  from  other  layers  so  that  we  can  drop  the  label  g.  In  practice 
much  the  same  result  is  found  for  aromatic  hydrocarbon  crystals,  where  the  bulk 
properties  are  dominated  by  thos*.  of  a  single  layer  once  allowance  is  made  for  the 
slab-shaped  crystal  implied  by  the  planewise  summation.* 

4.  DIELECTRIC  RESPONSE 

Given  the  foregoing  conclusions,  the  dielectric  response  of  a  film  is  adequately 
represented  by  that  of  a  single  layer.  The  first  quantity  of  interest  is  the  local  field 
F.  This  is  related  to  the  applied  field  by  the  dipole  tensor  T,  but  to  the  more  useful 
macroscopic  field  E  by  the  Lorentz-factor  tensor  L  =  T  -1-  «n,  where  n  is  the  unit 
vector  normal  to  the  layers."  Then  in  the  present  case  except  that  = 

7L  -f  1,  and  hence  Tr  L  =  1.  We  define  the  local-field  tensor  by  F  =  d  •  JS,  and 
obtain 


d  =  (1  -  L  •  a)-',  (1) 

where  a  =  alt(,v  is  a  dimensionless  reduced  polarizability.  The  calculated  local-field 
components  for  the  reduced  bead  polarizability  a/toV  =  0.2  are  given  in  Table  II. 
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1.109 
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1.109 
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1.018 

20° 

1.112 

0 

0.025 
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1.100 

0 

0.001 

0 

1.085 

40° 

1.085 

0 

-0.142 

0 

1.089 

0 

0.028 

0 

1.151 

Table  II.  Local-field  tensor  components  as  a  function  of  tilt  expressed  relative  to 
the  cartesian  axes  x  y  z  defined  in  Table  I. 

The  local  fields  are  not  greatly  different  from  the  macroscopic  field.  The 
dependence  on  tilt  is  not  very  marked,  and  is  not  monotonic  because  tilt  affects 
corresponding  components  of  L  and  a  in  opposite  senses  in  the  xz  plane. 


250  Model  Langmuir-Blodgett  Films 


The  linear  susceptibility  follows  from  d  as 

jr"*’  =  a  •  d  (2) 

and  the  relative  permittivity  is  then  e  =  jf**’  +  1 .  The  refractive  indices  are  obtained 
from  the  indicatrix  or  inverse  relative  permittivity  for  a  direction  e,  the  refractive 
index  n,  is  given  by  Hn^  =  e-€~'-e.  Calculated  refractive  indices  and  the  angle  </> 
between  the  principal  axes  of  refracdon  and  the  x  and  z  axes  are  given  in  Table  III. 


e 

n. 

4> 

0° 

1.105 

1.105 

1.421 

0° 

20° 

1.134 

1.096 

1.389 

9.3° 

o 

O 

L._ 

1.214 

1.104 

1.281 

31.2° 

Table  HI.  Refractive  indices  in  the  x,  y  and  z  direcdons  and  the  angle  between 
principal  axis  of  refraction  and  x  direction  as  a  funcdon  of  dlt. 


As  the  molecular  long  axis  tilts  away  from  the  vertical  z  axis  and  towards  the 
-X  axis,  /ij  decreases  and  increases  as  expected,  while  n,  increases  slightly.  The 
angle  <^>,  measured  from  +x  towards  +z,  increases  with  the  tilt  but  lags  about  10° 
behind  it,  so  that  the  principal  axis  does  not  exactly  follow  the  molecular  long  axis. 
The  quadradc  suscepdbility  is  given  by 

•  yff  :  dd  /  6oV  (3) 


where  the  superscript  T  denotes  the  transpose.  Calculated  components  of  are 
given  in  Table  IV. 
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Table  TV.  Quadradc  susceptibility  components  10“’  esu  as  a  function  of  tilt;  only 
independent  nonzero  components  are  shown. 


The  calculated  behaviour  of  Xin  be  understood  in  terms  of  the  negadve 


R.  W.  Munn  and  M.  M.  Shabat  251 


sign  of  the  ZZZ  component  of  /?  and  the  positive  sign  of  the  other  components,  but 
generally  the  smaller  off-diagonal  elements  of  P  enhance  the  effect  of  ^2zz- 
induces  additional  nonzero  components  of  so  that  at  40°  tilt  the  pattern  is  very 
different  from  that  at  zero  tilt,  with  no  longer  the  largest  component  in  magnitude 
or  even  the  most  negative  component.  Clearly  this  would  have  a  significant  effect  on 
nonlinear  optical  measurements. 


5.  DISCUSSION 

The  present  results  show  how  one  can  proceed  from  suitable  models  of  LB 
film  structure  and  of  molecular  response  to  calculate  the  linear  and  nonlinear  optical 
properties  of  the  film.  The  values  presented  here  are  illustrative  rather  than  definitive, 
but  they  do  yield  plausible  results.  The  models  are  simplified,  particularly  for  the 
distribution  of  molecular  response,  which  is  unlikely  to  be  the  same  for  each  "bead" 
in  molecules  of  practical  interest,  but  more  realistic  treatments  can  be  carried  out 
when  required.  As  indicated  in  the  Introduction,  these  calculations  can  be  combined 
with  molecular  dynamics  modelling  to  allow  one  to  predict  the  response  of  films 
made  from  specifically  designed  molecules,  and  hence  to  implement  a  systematic 
approach  to  this  area  of  molecular  electronics. 
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1  INTRODUCTION 

Organic  materials  are  attractive  for  nonlinear  optics  because  they  offer  the 
possibility  of  tailoring  the  molecules  to  modify  their  properties.  However,  this 
implies  but  also  that  one  understands  how  the  desired  material  properties  arise 
from  the  molecular  properties.  The  basic  understanding  is  well  established:  the 
material  response  is  given  by  the  molecular  response  transformed  from  molecular 
to  material  axes,  modified  by  local-field  factors  that  relate  the  polarizing  field  in 
the  material  to  the  macroscopic  electric  field  at  the  relevant  frequencies.  ^  Often 
the  local-field  factors  are  approximated  by  the  Lorentz  local-field  factor  using  the 
mean  refractive  index,  in  the  spirit  of  the  oriented-gas  model.^  In  crystals,  this 
approximation  is  known  to  be  reasonable  for  compact  molecules  but  not  for 
elongated  ones,^  and  hence  care  needs  to  be  taken  in  interpreting  and  predicting 
NLO  response  of  ordered  materials  composed  of  markedly  anisotropic  molecules. 

These  reservations  clearly  apply  to  Langmuir-Blodgett  films,  formed  from 
amphiphilic  molecules  with  a  hydrophilic  "head"  and  a  long  hydrophobic  "tail". 
The  effect  of  the  elongated  structure  on  the  local  field  is  not  obvious,  particularly 
if  the  molecular  axis  is  tilted  away  from  the  normal  to  the  film,  as  often  seems 
to  occur. Previous  algebraic  results  show  how  a  planewise  approach  can  be  used 
to  calculate  the  linear  and  nonlinear  optical  response  of  layered  materials  such 
as  LB  films.^  Here  we  report  numerical  calculations  for  simplified  model  LB 
films  in  order  to  explore  how  the  response  varies  with  molecular  tilt  and  lattice 
distortion.  This  work  forms  pan  of  a  coordinated  programme  to  design,  prepare 
and  characterize  LB  films  with  high  NLO  activity.® 


2  MODEL  STRUCTURE 

The  model  is  composed  of  molecules  treated  as  a  string  of  identical  "beads"  to 
represent  their  size,  shape  and  orientation.  The  parent  strucmre  has  the  molecules 
in  each  layer  of  the  film  packed  in  a  regular  triangular  array  with  their  long  axes 
perpendicular  to  the  layer.  Centring  alternate  layers  so  that  molecules  lie  over  the 
triangle  formed  by  molecules  in  the  preceding  layer  produces  a  hexagonal 


structure.  A  plan  view  of  the  structure  is  shown  in  Figure  1. 


Figure  1.  Plan  view  of  layer  showing  bases  of  the 
hexagonal  (light)  and  monoclinic  (bold)  cells. 


Apart  from  the  conventional  hexagonal  unit  cell,  the  structure  can  also  be 
described  by  the  centred  monoclinic  cell  indicated.  Its  short  edge  length  is  the 
nearest-neighbour  (NN)  distance  D,  and  its  long  edge  length  is  the  next-nearest- 
neighbour  (NNN)  distance  a/Jd,  where  D  is  the  diameter  of  a  "bead".  When 
molecules  tilt,  it  is  necessary  to  specify  not  only  the  tilt  but  also  the  plane  of  tilt. 
It  appears  that  NNN  tilt  prevails,  and  the  monoclinic  cell  is  then  preferable. 
There  are  also  indications  from  atomic-force  microscopy  that  some  LB  films  have 
a  distorted  hexagonal  structure,^  which  is  conveniently  modelled  by  changing  the 
axial  ratio  of  the  monoclinic  cell.  For  multilayers,  it  is  necessary  to  specify  the 
relationship  between  tilts  in  successive  layers,  which  may  be  the  same  or 
alternating,  as  shown  in  Figure  2. 


Figure  2.  Elevation  of  structures  showing  successive  layers  without  tilt,  with  the 
same  tilt,  and  with  alternating  tilt. 


3  METHOD 

Interactions  between  molecules  are  expressed  in  terms  of  planewise  dipole  sums 
giving  the  field  at  a  molecule  in  one  layer  due  to  dipoles  on  the  molecules  in  that 
or  another  layer.  These  sums  have  been  calculated  using  standard  expressions.*’^ 
They  are  calculated  separately  between  all  pairs  of  "beads"  and  then  averaged  to 
give  molecular  sums  T^(g),  where  a  and  0  are  cartesian  components  and  g  is 
the  number  of  layer  spacings  between  the  origin  molecule  and  the  layer  of 
dipoles.  Hence  g  =  0  denotes  interactions  within  a  layer,  g  =  1  interactions 
between  adjacent  layers,  and  so  on.  The  sums  are  scaled  by  a  factor  v/4t,  where 
V  is  the  molecular  volume,  to  give  dimensionless  quantities. 


As  reported  below,  interactions  for  ^  >  1  are  very  small,  so  that  it 
suffices  to  consider  only  the  sum  T(0)  and  the  corresponding  planewise  Loren tz- 
factor  tensor  L  =  T(0)  +  /in,  where  n  is  the  normal  to  the  layers.  The  local 
polarizing  electric  field  F  is  related  to  the  macroscopic  field  E  by 

F=(1  -L.ari*£  =  d-£,  (1) 

where  a  =  ale^v  is  a  dimensionless  reduced  polarizability  and  d  is  the  local-field 
tensor.  The  relative  permittivity  is 

e  =  1  +  a-d  =  1  +  (2) 

where  is  the  linear  susceptibility,  and  in  direction  e  the  refractive  index 
is  given  by  1/n/  =  e’€~^’e. 

The  quadratic  susceptibility  is^’^® 

/2)  =  dT.b:dd,  (3) 

where  the  superscript  T  denotes  the  transpose  and  b  =  is  a  reduced 
molecular  first  hyperpolarizability.  The  cubic  susceptibility  comprises  a  direct 
term  and  a  cascading  term 

y3)  _  d’^.c-ddd  +  2d^*b:dd*d*L*b:dd  s  +  x^^^\  (4) 

where  c  =  y/^QV  is  a  reduced  molecular  second  hyperpolarizability. 

The  molecular  response  is  treated  in  a  simplified  way  designed  to  make 
changes  of  film  structure  the  dom.inant  influence  on  film  response.  In  molecular 
axes  .Y  y  Z,  the  polarizability  is  taken  as  {a,  a,  sa),  with  a/egV  fixed  as  0.2.  The 
first  hyperpolarizability  is  taken  to  have  nonzero  components  diagonal  in 
AB,  with  components  (jS,  0,  -  s0)  diagonal  in  AB  and  jS/egV  fixed  as  3.8  pm  V'^ 
The  second  hyperpolarizability  is  taken  to  have  only  the  component  =  sy 
nonzero,  with  y/eQV  fixed  as  200  (pm  V'^)^.  These  choices  are  guided  by 
published  calculations  of  molecular  response. 


4  RESULTS 

Planewise  dipole  tensor  sum  components  for  g  =  0  and  1  as  a  function 
of  tilt  are  given  in  Table  1  for  j  =  5  relative  to  x  y  z  cartesian  axes,  where  z  is 
normal  to  the  layers  and  xz  is  the  plane  of  tilt.  For  g  =  1,  the  values  are  at  most 
1  %  of  the  corresponding  components  for  g  =  0,  and  hence  terms  for  g  >  1  can 
be  neglected.  For  g  =  0  at  zero  tilt,  the  values  give  a  diagonal  Lorentz-factor 
tensor  close  to  (0.5,  0.5,  0),  the  depolarization  factor  expected  for  a  needle- 
shaped  Lorentz  cavity,  consistent  with  the  elongated  molecular  shape.  Tilt  affects 
little,  but  brings  and  L„  close.  Very  similar  results  are  obtained  for 


Table  1.  Planewise  sum  components  T^^{g)  for  various  tilts  6. 


e 

a0 

XX 

xz 

>7 

zz 

g  =  0 

0° 

0.491 

0 

0.491 

-0.982 

20° 

0.440 

0.093 

0.467 

-0.907 

o 

O 

0.281 

0.148 

0.420 

-0.701 

g  =  i 

0° 

-0.0016 

0 

-0.0016 

0.0031 

o 

O 

-0.0010 

-0.0023 

-0.0007 

-0.0002 

o 

O 

0.0016 

0.0011 

0.0043 

-0.0059 

values  of  s  up  to  10,  and  so  future  results  all  refer  to  j  =  5.  The  effect  of 
distortion  was  explored  by  expanding  one  edge  of  the  monoclinic  cell  base  and 
contracting  the  other  in  the  same  proportion.  Distortions  up  to  10%  produced 
linear  variations  in  the  planewise  sum  components  of  similar  proportions. 

The  local  field  component  depends  on  tilt  only  through  L  and  so 
varies  little,  from  1. 109  at  zero  tilt  to  1.092  at  40°  tilt.  In  the  x  and  z  directions, 
tilt  affects  the  Lorentz-factor  tensor  and  polarizability  components  in  opposite 
senses,  so  there  is  also  no  great  change.  The  local  field  component  increases 
from  1.109  at  zero  tilt  to  1.118  at  20°  but  then  decreases  to  1.100  at  40°,  while 
<i_  increases  from  1.018  at  zero  tilt  to  1.064  at  20°  and  1.166  at  40°.  Off- 

ZZ 

diagonal  components  are  always  less  than  0.05. 

The  refractive  index  ity  varies  little  with  tilt  because  the  polarizability 
along  y  does  not  change.  The  variation  of  the  other  refractive  indices  largely 
reflects  that  of  the  polarizability:  increases  from  1.105  at  zero  tilt  to  1.217  at 
40°,  and  n,  decreases  from  1.421  to  1.285.  The  principal  directions  in  the  xz 
plane  remain' within  1°  of  the  molecular  axes.  Distortion  of  the  lattice  as 
described  earlier  has  little  effect  on  the  refractive  indices:  10%  distortion  at  zero 
tilt  changes  by  -h0.0009,  riy  by  -0.0007,  and  n.  by  -0.0018. 

Quadratic  susceptibility  components  as  a  function  of  tilt  are  given  in  Table 

2.  The  assumed  form  of  0  ensures  Kleinman  symmetry.  The  effects  of  tilt  are 
marked,  and  monotonic  except  for  Xxxx  where  changes  sign  as  6  increases. 

The  direct  parts  of  the  cubic  susceptibility  components  are  given  in  Table 

3.  Tilt  decreases  the  zzzz  component  and  increases  the  others  in  magnitude. 
Cascading  parts  are  also  given  in  Table  3  with  P  treated  as  having  only  an  axial 


Table  2.  Quadratic  susceptibility  components  for  various  tilts  0. 


9 

a0y 

XXX 

xxz 

xyy 

xzz 

zyy 

zzz 

0° 

0 

4.9 

0 

0 

4.9 

-20.9 

20° 

-3.8 

0.2 

-1.8 

10.7 

4.8 

-18.5 

o 

O 

0.3 

-8.8 

-3.2 

13.5 

■H 

-8.3 

component,  like  y-  These  parts  reinforce  the  direct  ones,  except  for  x-t-  where 
they  tend  to  offset  the  decreases  with  tilt,  and  are  typically  a  quarter  of  the 
magnitude  at  20°  tilt  and  over  half  at  40°.  Including  all  components  of  /? 
complicates  the  pattern  by  introducing  numerous  additional  components,  the 
largest  being  Xyyyy  which  reaches  5700  {pm  at  40°  tilt,  and  modifying 
others,  notably  x^^  which  reaches  11  900  (pm  V^)^. 

Table  3.  Cubic  susceptibility  components  Xa/j^^/Cpm  V'^)^  for  various  tilts  9. 
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xxxx 

xxxz 

xxzz 

xzzz 

zzzz 

Direct  contribution 


0° 

0 

0 

0 

0 

10800 

20° 

200 

-600 

1500 

-3900 

10500 

o 

O 

2100 

-2800 

3600 

-4800 

6300 

Cascading  contribution 


0° 

0 

0 

0 

0 

200 

20° 

100 

-200 

400 

-1100 

2900 

o 

O 

1300 

-1700 

2300 

-3000 

4000 

5  DISCUSSION 

The  numerical  values  used  in  the  present  work  are  arbitrary  but  plausible.  Hence 
a  number  of  the  qualitative  results  should  be  fairly  widely  applicable. 

For  elongated  molecules,  dipolar  interactions  between  different  layers  are 
essentially  negligible,  even  for  a  tilt  as  large  as  40°.  This  greatly  simplifies  the 
dielectric  theory.  The  local  fields  are  surprisingly  little  affected  by  tilt,  so  that 
the  linear  susceptibility  and  refractive  indices  are  determined  largely  by  the 


relevant  components  of  the  effective  polarizability.  Birefringence  measurements 
could  thus  provide  reliable  information  on  tilt,  but  probably  not  on  lattice 
distortion. 

The  nonlinear  susceptibilities  can  be  fairly  readily  interpreted  in  terms  of 
the  effect  of  tilt  on  the  molecular  polarizability  referred  to  the  film  x  y  z  axis 
system.  The  lower  symmetry  of  the  tilted  structures  produces  a  much  more 
complex  set  of  components  and  correspondingly  requires  a  more  complicated 
interpretation  of  measurements,  especially  since  components  induced  by  tilt  can 
be  comparable  with  those  for  zero  tilt.  In  the  cubic  susceptibility,  the  cascading 
terms  can  be  as  large  as  the  direct  terms,  so  that  a  simple  interpretation  based 
solely  on  the  second  hypeipolarizability  may  be  misleading. 

Future  work  is  required  on  more  realistic  molecular  structures  in  which 
the  response  is  not  assumed  to  be  uniformly  distributed  but  is  concentrated  in  the 
NLO  chromophore. 
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